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INTRODUCTION 



The fmiction Q 



In this paper we use the initial-state-radiation (ISR) 
technique to study the e^e~ — >■ pp process in a wide range 
of center-of-mass (cm.) energies. The Born cross section 
for the ISR process e+e^ pp^ (Fig. [U, integrated over 
the nucleon momenta, is given by 



d Ce+e- — J-PP7 (-^pp) 

dMpp d cos 9* 



9 /If 

^^W{s,x,9*)app{Mpp), (1) 



where <Jpp{m) is the Born cross section for the nonradia- 
tive process e^e~ — pp, Mpp is the pp invariant mass, 
is the nominal e~^e~ cm. energy, x = 2E* 



1 — Mpp/s, and E* and 9* are the ISR photon energy 
and polar angle, respectively, in the e"^e~ cm. frame. ^ 
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W{s,x,9;) = ^l^-^ 
TTx \ sm ( 



(2) 
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is the probability of ISR photon emission for 9* ^ 
rrie/y/s, where a is the fine-structure constant and me is 
the electron mass. The cross section for the e+e" — > pp 
process is given by 



3 Mpp 



\GMiMpp)\' 



2ml 

-f|GB(Mpp)p 
i^ipp 



where 



is the nominal proton mass, /3 



(3) 



pp, and C is the Coulomb correction fac- 



1 - 4m2/Af, 

tor (see, for example, Ref. [H and references therein), 
which makes the cross section nonzero at threshold (C = 
y/(l — e~y) with y = ira/ [3). The cross section depends 
on the magnetic form factor {Gm) and the electric form 
factor (Ge)', at threshold, \Ge\ = \Gm\- From the mea- 
surement of the cross section a linear combination of the 
squared form factors can be determined. We define the 
effective form factor 



I (Mpp) I 



'\Gm{MppW + 2ml/M^p\GE{MppW 



1 + 2ml/ M^p 



(4) 

which is proportional to the square root of the measured 
e'^e^ pp cross section. 



6 




FIG. 1: The Feynman diagram for the signal ISR process 
e~ VP')- 



The modulus of the ratio of the electric and magnetic 
form factors can be determined from the distribution of 
the angle between the proton momentum in the pp 
rest frame and the momentum of the pp system in the 
e^e~ cm. frame. This distribution can be expressed as 
a sum of terms proportional to |G7\/p and |G£;p. The full 
differential cross section for e^e^ ppj can be found, 
for example, in Ref. Q. The 6p dependences of the Ge 
and Gm terms are close to sin^ 6p and 1 + cos^ 9p, respec- 
tively. 

Direct measurements of the e"^e~ — > pp cross section 
have been performed in e+e~ experiments P-[l0j. Most 
of these measurements have an accuracy of (20-30)%. 
The cross section and the proton form factor were de- 
duced assuming jG^'j = |Gm|, and the measured proton 
angular distributions 0, did not contradict this 
assumption. More precise measurements of the proton 
form factor have been performed in — >■ e"'"e~ experi- 
ments [ni-El- In the PS170 experiment at LEAR, 
the proton form factor was measured from threshold {pp 
annihilation at rest) up to a cm. energy 2.05 GeV. The 
ratio \Ge/Gi\i\ was measured with about 30% accuracy 
and was found to be compatible with unity. The LEAR 
data show a strong dependence of the form factor on 
cm. energy near threshold, and very little dependence in 
the range 1.95-2.05 GeV. Analyses from Fermilab exper- 
iments E760 and E835 show a strong decrease 
of the form factor at cm. energies higher than 3 GeV, in 
agreement with perturbative QCD, which predicts that 
the dependence should be 0^(111?) /m'^ . However, the re- 
cent precision e~^e~ measurement based on CLEO 
data indicates that the decrease of the form factor at 
energies about 4 GeV is somewhat slower. 

The previous BaBar^ study [3] of the process e+e^ ^ 
pp using the ISR technique was based on about half of 
the data that were finally collected in the experiment. 
The e'^e~ — >■ pp cross section was measured for cm. en- 
ergies up to 4.5 GeV. This measurement yielded a signifi- 
cant improvement in precision for energies below 3 GeV. 



In contrast to previous e~^e~ and pp experiments, the 
BaBAR measurement did not assume that |G_b| = \Gm\- 
The ISR approach provides full 9p coverage, and hence 
high sensitivity to \Ge/Gm\- The energy dependence of 
the form-factor ratio \Ge/Gm\ was measured for cm. en- 
ergies below 3 GeV. For energies up to 2.1 GeV, this ratio 
was found to be significantly greater than unity, in dis- 
agreement with the PS 170 measurement [ll[. 

In this work we update the analysis of Ref. [ij] using 
the full BaBAR data sample collected at and near the 
T(4S) resonance. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

We analyse a data sample corresponding to an inte- 
grated luminosity of 469 fb^^ recorded with the BaBAR 
detector [l^ at the SLAG PEP-II asymmetric-energy col- 
lider. At PEP-II, 9-GeV electrons collide with 3.1-GeV 
positrons at a cm. energy of 10.58 GeV (the T(4S) mass). 
About 10% of the data are collected at 10.54 GeV. 

Charged-particle tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH), operating in a 1.5-T axial magnetic field. The 
transverse momentum resolution is 0.47% at 1 GeV/c. 
The position and energy of a photon-produced cluster 
are measured with a CsI(Tl) electromagnetic calorime- 
ter, which yields an energy resolution of 3% at 1 GeV. 
Charged-particle identification is provided by specific 
ionization (dE/dx) measurements in the SVT and DCH, 
and by an internally reflecting ring-imaging Cherenkov 
detector (DIRC). Muons are identified in the solenoid's 
instrumented fiux return (IFR), which consists of iron 
plates interleaved with either resistive plate chambers or 
streamer tubes [l6| . 

Signal and background ISR processes are simulated 
with Monte Carlo (MC) event generators based on 
Ref. [l3|- The differential cross section for e+e" — ppj 
is taken from Ref. To analyze the experimental pro- 
ton angular distribution, two samples of signal events 
are generated, one with Ge — and the other with 
Gm ~ 0. Since the polar-angle distribution of the ISR 
photon is peaked near 0° and 180°, the MC events are 
generated with a restriction on the photon polar an- 
gle: 20° < 9* < 160° (the corresponding angular range 
in the laboratory frame is 12° < 9j < 146°). Extra 
soft-photon radiation from the initial state is generated 
by the structure function method [l^- To restrict the 
maximum energy of the extra photons, the invariant 
mass of the hadron system and the ISR photon is re- 
quired to be greater than 8 GeV/c^. For background 
e'^e^ — > ^^/^~7, 7r"'"7r~7, and K~^K~^ processes, final- 
state bremsstrahlung is generated using the PHOTOS 
package [l9l |. Background from e^e~ qq is simulated 
with the JETSET HO] event generator; JETSET also 
generates ISR events with hadron invariant mass above 
2 GeV/c^, and therefore can be used to study ISR back- 
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ground with baryons in the final state. The dominant 
background process, e+e" — )■ pp7r°, is simulated sepa- 
rately by generating the angular and energy distributions 
for the final-state hadrons according to three-body phase 
space. 

The response of the BaBAR detector is simulated us- 
ing the Geant4 [2lj package. The simulation takes into 
account the variations in the detector and beam back- 
ground conditions over the running period of the experi- 
ment. 



III. EVENT SELECTION 

The preliminary selection of e+e^ — > pp7 candidates 
requires that all of the final-state particles be detected 
inside a fiducial volume. Since a significant fraction of 
the events contains beam-generated spurious tracks and 
photon candidates, we select events with at least two 
tracks with opposite charge and at least one photon can- 
didate with E* > 3 GeV. The polar angle of the photon 
is required to be in the well-understood region of the 
calorimeter: 20° < 9^ < 137.5°. Each charged-particle 
track must originate from the interaction region, have 
transverse momentum greater than 0.1 GeV/c and be 
in the angular region 25.8° < 9 < 137.5°. The latter 
requirement is needed to provide particle identification 
(PID) from the DIRC. To suppress background from ra- 
diative Bhabha events, we reject events for which the ra- 
tios of calorimetric energy deposition to momentum for 
the two highest-momentum tracks satisfy the condition 

(i^cau/pi - + (-Bcal,2/P2 " 1)' < 0.35^. 

For events passing the preliminary selection, a kine- 
matic fit is performed to the e^e^ — > h^h'^j hypoth- 
esis with requirements of total energy and momentum 
conservation. Here h can be e, ^, tt, K or p, and 7 
is the photon candidate with the highest energy in the 
e"*'e~ cm. frame. For events with more than two charged 
tracks, the fit uses the two oppositely charged tracks that 
pass closest to the interaction point. The MC simulation 
does not accurately reproduce the shape of the photon 
energy resolution function. This leads to a difference in 
the distributions of the oi the kinematic fit for data 
and simulated events. To reduce this difference, only 
the measured direction of the ISR photon is used in the 
fit; its energy is a free fit parameter. For each of the five 
mass hypotheses, the corrected angles and energies of the 
particles and the value are obtained from the fit. 

The expected number of events from the background 
processes e^e~ — > 7r+7r~7, /i"'"/i~7, and exceeds 
the number of signal events by two to three orders of mag- 
nitude. To suppress these backgrounds, we require that 
both charged particles be identified as protons accord- 
ing to the specific ionization (dE/dx) measured in the 
SVT and DCH, and the Cherenkov angle measured in 
the DIRC. This requirement suppresses pion and muon 
backgrounds by a factor of 3 x 10"*, and kaon background 
by a factor 10"*, with a loss of approximately 30% of the 
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FIG. 2: The Xp distribution for simulated e^e~ — > pfry (solid 
histogram) and e"*"e~ — > K^K^'y (dashed histogram, arbi- 
trary normalization) events. 



signal events. 

Background is further suppressed through require- 
ments on the x? of the kinematic fit: Xp < 30 and 
Xji > 30, where Xp and Xk the values of the kine- 
matic fit for the proton and kaon mass hypotheses, re- 
spectively. The Xp distribution for simulated ppj events 
is shown in Fig. [2l The long tail at high x^ is due to 
events with extra soft photons emitted in the initial state. 
The dashed histogram is the Xp distribution for K^K~j 
simulated events. The x^ requirements lead to the loss 
of 25% of the signal events, but provide additional back- 
ground suppression by a factor of 50 for pion and muon 
events, and a factor of 30 for kaon events. 

The pp invariant mass distribution is shown in Fig. [3] 
for the 8298 data events that satisfy the selection criteria. 
Most of the events have invariant mass below 3 GeV/c^. 
Signals from J ftp — > pp and ^^{23) — > pp decays are 
clearly seen. 



IV. BACKGROUND SUBTRACTION 

Potential sources of background in the sample of se- 
lected e+e~ — )• pp^ candidates are the processes e+e^ ^■ 
7r'*"7r~7, e"*"e~ — > K'^ K~^, e^e~ — >■ /i"'"^~7, and e+e~ — >■ 
e+e~7 in which the charged particles are misidentified as 
protons. Background contributions from processes with 
protons and neutral particle(s) in the final state, such as 
e+e^ — > ppTT^, ppr/, ppi:^^, etc., are also anticipated. 

Of particular interest is the possible background from 
the process e+e~ — >■ pp"f with the photon emitted from 
the final state. Due to the different charge-conjugation 
parity of the amplitudes corresponding to initial-state ra- 
diation and final-state radiation (FSR) , their interference 
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FIG. 3: The pp invariant mass spectrum for the selected data 
PP7 candidates. The left edge of the plot corresponds to the 
pp threshold. 

does not contribute to the total e+e" — ppj cross sec- 
tion. The contribution of the FSR amplitude is estimated 
to be [m da/dm « \ F^y,\^8ma^ 13/ {27 s^), where Fax is 
the axial proton form factor. Assuming |-Fax| ~ |Gm|, the 
ratio of the FSR and ISR cross sections is estimated to 
be about 10"'^ for pp invariant masses below 4.5 GeV/c^. 
We conclude that the FSR background is small and so 
may be neglected. 

A. Background contributions from 

and e~^e~ e~^e~y 

To estimate the background contribution from e~^e~ — >■ 
7r"*'7r~7, data and simulated tt'^'!t~j events are selected 
with the following requirements on PID and on the of 
the kinematic fit: 

1. one proton candidate, < 20; 

2. one proton candidate, < 30, Xk > 30; 

3. two proton candidates, xi < 20; 

4. two proton candidates, Xp < 30, Xk 30. 

Here xi is the x^ of the kinematic fit for the pion mass 
hypothesis. The fourth set of conditions corresponds to 
the standard selection criteria for ppj candidates. 

The invariant mass Mtttt of the two charged particles 
under the pion-mass hypothesis is calculated; the M^^ 
distributions for data selected with criteria 2 and 3 are 
shown in Fig. 21 The data spectra are fit with a sum of 
the mass spectra for simulated 7r"''7r~7 events (p-meson 
line shape with w-p interference) and a linear background 
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FIG. 4: (a) The M^r spectrum for data events with Xp < 30 
and Xk > 30, and one proton candidate (selection 2 in the 
text); (b) the same spectrum for data events with < 20 
and two proton candidates (selection 3 in the text). The 
histograms are the results of the fit described in the text. 



term. The numbers of tttt"/ events with 0.5 < A/tttt < 1 
GeV/c^ obtained from the fits for selections 1-3 are listed 
in TablelH together with the corresponding numbers from 
the 7r+7r~7 MC simulation. The spectrum for selection 
4 is fit by a linear function; no p-mcson contribution is 
needed to describe this spectrum. 

Since the simulation correctly predicts the numbers of 
pion events for selections 1-3, it can be used to estimate 
the pion background for the standard selection 4. We 
observe no events satisfying the standard selection cri- 
teria in the 7r7r7 MC sample. The corresponding upper 
limit on ttttj background in the data sample is 5.2 events 
at 90% confidence level (CL). The estimated pion back- 
ground is less than 0.1% of the number of selected ppj 
candidates. 
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TABLE I: The numbers of nnj events for data and MC simu- 
lation with 0.5 < A'l^TT < 1 GeV/c^ that satisfy different selec- 
tion criteria for data and MC simulation. The data numbers 
are obtained from the fits to the M^^ distributions described 
in the text. 



selection 


data 


MC 


f 


15310 ± 160 


14800 ± 180 


2 


400 ± 60 


460 ± 30 


3 


41 ±8 


48± 11 



Similarly, the number of e'^e^ — > K'^ K~j events can 
be estimated from the number of events in the </> meson 
peak in the distribution of invariant mass of the charged 
particles calculated under the kaon hypothesis. It is 
found that the K^K^j MC simulation predicts reason- 
ably well the numbers of kaon events in the data sample 
with one identified kaon and the standard conditions, 
and in the data sample with two identified kaons and 
Xk < 20. Therefore we use the MC simulation to esti- 
mate kaon background for the standard selection. The 
estimated background, 1.6 ± 0.8 events, is significantly 
less than 0.1% of the number of data events selected. 

The kinematic properties of the e+e^ — > e+e~7 pro- 
cess are used to estimate the electron background. About 
50% of e^e~7 events have e"*"e~ invariant mass between 3 
and 7 GcV/c^ and cos?/)* < —0.98, where -0* is the angle 
between the two tracks in the initial e"^e~ cm. frame. In 
the event sample with two proton candidates we do not 
find events having the above characteristics. The corre- 
sponding 90% CL upper limit on the e"'"e~7 background 
in the data sample is 4.6 events (2 events with Mpp < 4.5 
GeV/c2). 

To compare MC simulation and data for the process 
e~^e~ — >■ fi^ fi~j, we use a subsample of events selected 
with the requirement that both charged particles be iden- 
tified as muons. Muon identification is based on IFR 
information, and does not use DIRC or dE/dx informa- 
tion, which are necessary for proton identification. In 
the data samples with one or two identified protons ob- 
tained with the standard selection, we select 86 and 2 
muon-identified events, respectively. These numbers can 
be compared with 60 ± 16 and zero events expected from 
the e+e~ — > /i+/^~7 simulation. Taking into account that 
the ratio of the total number of /x+/x~7 events to those 
with two identified muons is about two-to-one, we esti- 
mate the /Lt^/z~7 background for the standard selection 
criteria to be 4.0 ± 2.8 events. 

The combined background from the processes e+e^ 
C~^C~'y, C = TT, K,e, ^ is less than 0.2% of the number 
of selected pp7 candidates, and so can be neglected. 



B. Background from e^e — ^ ppn'^ 

The dominant source of background to the e+e" -> 
PP7 process arises from e~^e~ — > ppir^ . A significant frac- 
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FIG. 5: (a) The Mpp spectrum and (b) the cos 9p distribution 
for selected e"'"e~ — >■ ppTv" candidates in data. In each figure, 
the shaded histogram shows the background contribution es- 
timated from the M-y-y sidebands. 



tion of ppTr" events with an undetected low-energy pho- 
ton, or with merged photons from the tt" decay, is recon- 
structed under the ppj hypothesis with a low value of x^, 
and so cannot be separated from the process under study. 
This background is studied by selecting a special subsam- 
ple of events containing two charged particles identified 
as protons and at least two photons with energy greater 
than 0.1 GcV, one of which must have cm. energy above 
3 GeV. The two-photon invariant mass is required to be 
in the range 0.07-0.20 GeV/c^, which is centered on the 
nominal tt^ mass. A kinematic fit to the e'^e" — > ppj'f 
hypothesis is then performed. Requirements on the 
of the kinematic fit (x^ < 25) and the two-photon in- 
variant mass (0.1025 < M^^ < 0.1675 GeV/c^) are im- 
posed in order to select e+e~ — >■ ppir^ candidates. The 
sidebands 0.0700 < M^^ < 0.1025 GeV/c^ and 
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FIG. 6: The cos6^ distribution for e^e~ — >■ pp7r° event can- 
didates for data (points with error bars) and simulation (his- 
togram) . 

0.1675 < M^^ < 0.2000 GeV/c^ are used to estimate 
background. The Mpp spectra and cos^p distributions 
for data events from the signal and sideband M^j regions 
are shown in Fig. [SJ The total number of selected events 
is 148 in the signal region and 12 in the sidebands. The 
number of e+e^ ppn^ events in the M^^-y sidebands 
expected from MC simulation is 5.4. 

The ppTT^ selection criteria described above are applied 
to simulated e+e^ — > qq events generated with the JET- 
SET program. In the mass region Mpp < 5 GeV/c^, 
the predicted number of e+e" — !• ppir'^ events is 120 ± 9. 
These events exhibit an enhancement in the Mpp distri- 
bution near the pp threshold, similar to that in data. 
However, the cos^p distribution of MC events is peaked 
near cos 9p = ±1, whereas for the data the distribution is 
flat. To study the e+e" — >■ ppir'^ background, the sample 
of simulated e+e^ ppn'^ events is generated according 
to three-body phase space, but with an additional weight 
proportional to {Mpp — 2mp)^/^ to imitate the AIpp distri- 
bution observed in data. The resulting generated cos 9p 
distribution is flat. 

In Fig. |6] the cos 9^ distribution for selected data and 
simulated e"*'e~ — > ppw'^ events is shown, where 0* is 
the tt" polar angle in the e~^e~ cm. frame. It is seen 
that the data and simulated distributions differ slightly. 
Since we do not observe a significant variation of the 
cos 0* distribution with Mpp in data, we use the data 
distribution averaged over Mpp (Fig. [6]) to reweight the 
e"'"e~ — > ppir'^ simulation. 

From the reweighted simulation, we calculate the ratio 
(-?^Mc) of the Mpp distribution for events selected with 
the standard criteria to that selected with the ppn^ 
criteria. The value of the ratio KyiciMpp) varies from 
3.4 near the Mpp threshold to 2.0 at 5 GeV/c^. The ex- 
pected Mpp spectrum for the e~^e~ ppir'^ background 



FIG. 7: The expected Mpp spectrum for e^e~ — > pp7r° events 
selected with the standard pfry criteria. The spectrum is ob- 
tained by scaling the data distribution shown in Fig. Ela) by 
the factor KMc{Mpp) described in the text. 

satisfying the pp^ selection criteria is shown in Fig. [7J 
and is evaluated as ifMc(-^pp) x {dN/dMpp)data, where 
{dN/dMpp)data is the mass distribution for e+e" ppn^ 
events obtained above (Fig. El^a)). The number of se- 
lected e+e~ ppj candidates and the expected num- 
ber of e+e^ ppi:^ background events are given for 
different pp mass ranges in Table [III The background 
contribution grows from 5% near pp threshold to 50% 
at Mpp K, 4 GeV/c^. All observed pp^ candidates with 
Mpp > 4.5 GeV/c^ are consistent with ppir'^ background. 

C. Other sources of background 

Other possible background sources are ISR processes 
with higher final-state multiplicity {e'^e^ — > ppir^j, 
pp27r°7, ...), and direct e'^e~ annihilation processes 
other than e'^e^ — >■ ppn'^ (e^e^ — >■ pprj, e^e^ — > pp2TT^, 
and so on). All of these processes arc simulated by JET- 
SET. The simulation leads to the prediction that the ISR 
background is 55 ± 6 events, and that the direct annihi- 
lation background is 40 ± 5 events. The total predicted 
background from these two sources is about 1.2% of the 
number of selected pp"f candidates. We do not perform a 
detailed study of these background processes. Their con- 
tribution is estimated from data by using the sideband 
region, as described below in Sec. IIVDI 

D. Background subtraction 

Table IIIII summarizes the expected number of back- 
ground events estimated in the previous sections. The 
"other ISR" and "e+e~" columns show the background 
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TABLE II: The number of selected pp7 candidates, Npp^, and the number of background events from the e^e — >■ ppn 
process, Npp^o, for different ranges of Mpp. The pp mass ranges near the J/tp and ?/'(25') resonances are excluded. 

Mpp (GeV/c^) < 2.50 2.50-3.05 3.15-3.60 3.75-4.50 > 4.5 

Npp-, 6695 592 76 29 9 

N^^^o 321 ±37 66 ± 15 26 ± 9 17 ±6 6±3 



TABLE III: The number of selected pp7 candidates from the mass region Mpp < 4.5 GeV/c^ with xl < 30 (iVi) and 30 < Xp < 60 
{N2) for signal and for different background processes; Pi is the ratio N2/N1 obtained from simulation. The first column shows 
the numbers of pp7 candidates selected in data. The numbers for e^e~ pp^ are obtained from data using the background 
subtraction procedure described in the text. 





data 


- 
ppn 


e+e- 


Other ISR 


PPl 


iVi 


8298 


448 ± 42 


40 ±5 


55 ±6 


7741 ±113 


iV2 


560 


79 ±7 


76 ±7 


74 ±7 


337 ± 16 


ft 




0.175 ±0.04 


1.88 ±0.29 


1.34 ± 0.18 


0.0435 ± 0.0020 



contributions estimated with JETSET that result from 
ISR processes, and from e'^e" annihilation processes 
other than e+e" — )■ ppn^. Because JETSET has not 
been verified precisely for the rare processes contributing 
to the candidate sample, we use a method of back- 
ground estimation that is based on the difference in 
distributions for signal and background events. The first 
and second rows in Table lllll show the expected numbers 
of signal and background events with Xp < 30 (A^i) and 
30 < Xp < 60 (iV2)- The last row lists /3i, the ratio of N2 
to Ni . From Table IIIH it is evident that the coefficient 
/3i for signal events, and for background events from the 
processes with higher hadron multiplicity ("e'''e~" and 
"Other ISR" columns), are very different. This differ- 
ence can be used to estimate the background from these 
two sources, as follows. First, the ppir'^ background deter- 
mined as described in Sec. IIVBI is subtracted from data. 
Then, from the corrected numbers of events in the sig- 
nal and sideband regions, N[ and N2, the numbers of 
signal and background (from "e"'"e~" and "ISR" sources) 
events with Xp < 30 can be calculated: 



N ■ 

Nhkg 



1 l^pp'y / ^bkg 



(5) 



where /3bkg is the ratio of the fractions of events in the 
sideband and signal x^ regions averaged over all back- 
ground processes of the "e+e~" and "ISR" types. For 
this coefBcient, f3bkg = 1.6 ± 0.3 is used; it is the average 
of /3e+e- and /3isr with the uncertainty (/?e+e- — /3isr)/2. 

In Table IIIIl it is also evident that ppj events domi- 
nate the x^ sideband region. Therefore, the background 
is very sensitive to the accuracy of the Ppp-y coefficient. In 
particular, the data-Monte Carlo difference in the dis- 
tribution can lead to a systematic shift of the result. The 
simulation of the distribution for ppj events is vali- 



dated using data and simulated e"' 



/i"*"/! 7 events. 



study, and can be selected with negligible background. 
The ratio of the /? coefficients for e+e^ /.i^n^-f data 
and simulation is independent of the /i^/i^ mass and is 
equal to 1.008 ± 0.008. This ratio is used to correct the 
Pppj value obtained from simulation, which varies from 
0.043 at pp threshold to 0.048 at 4.5 GcV/c^. 

With the method described above, the total numbers of 
e"'"e~ — > ppj events (A^sig) and background events from 
"e+e^" and "ISR" sources (Nbkg) in the signal region are 
found to be 7741±95±62 and 109±16±25, respectively. 
The main source of the systematic uncertainty on Ngig 
is the uncertainty in the ppn^ background. The number 
of background events is in good agreement with the esti- 
mate from simulation, (40 ± 5) ± (55 ± 6) = 95 ± 8. The 
total background in the Xp < 30 region is 531±51 events, 
which is about 7% of the number of signal events. 

The background subtraction procedure is performed in 
each pp mass interval. The number of selected events for 
each interval after background subtraction and unfolding 
event migration between intervals (see Sec. IVIip is listed 
in Table IVTl The events from J/tjj and ip{2S) decays are 
subtracted from the contents of the corresponding inter- 
vals (see Sec. lVin)) . 



V. ANGULAR DISTRIBUTIONS 

The modulus of the ratio of the electric and magnetic 
form factors can be extracted from an analysis of the dis- 
tribution of 9p, the angle between the proton momentum 
in the pp rest frame and the momentum of the pp system 
in the e'^e^ cm. frame. This distribution is given by 



dN 



dcos Or, 



These are very similar kinematically to the process under 



A Hm {cos Op, Mpp) 



He (cos Op, Mpp) . (6) 
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TABLE IV: The number of selected pp7 candidates (N) and 
the number of background events {Ntkg) for each pp mass 
interval; \Ge/Gm\ is the fitted ratio of form factors. 



are shown in Fig. [U the values of \Ge/Gm\ obtained are 
listed in Table IIVI and shown in Fig. [TUl The curve in 



Mpp, GeV/c^ 


iV 


Nbkg 


1.877-1.950 


1162 


19 ± 10 


1.950-2.025 


1290 


53 ±16 


2.025-2.100 


1328 


63± 14 


2.100-2.200 


1444 


118 ±28 


2.200-2.400 


1160 


126 ± 26 


2.400-3.000 


879 


122 ± 22 



\Ge/G 



M 



+0.16+0.06 
-0.14-0.05 



1 ofj+0. 15+0.05 

1.4i 

1 oQ+0. 15+0.07 
-'-■''•'-0.14-0.07 
1 9f;+0.14+0.10 
^■^"-0.13-0.09 
1 r,^+0.16+0.10 
^■"^-0.16-0.10 

1.04: 



-0.24- 
-0.25- 



0.15 
0.15 



The functions H m {cos 6p, Mpp) and HE{cosdp, Mpp) do 
not have an analytic form, and so arc determined using 
MC simulation. To do this, two samples of e+e~ — ?> pp7 
events are generated, one with Ge = and the other 
with Gm = 0. The functions obtained are close to the 
1 + cos^ 9p and sin^ dp functions describing angular dis- 
tributions for the magnetic and electric form factors in 
the case of e+e~ — )■ pp. 

The observed angular distributions are fit in six ranges 
of pp invariant mass from threshold to 3 GeV/c^. The fit 
intervals, the corresponding numbers of selected events, 
and the estimated numbers of background events arc 
listed in Table IIVI For each pp mass interval and each 
angular interval the background is subtracted using the 
procedure described in Section IIVDI The angular dis- 
tributions obtained are shown in Fig. [51 The distribu- 
tions are fit to Eq. ^ with two free parameters: A (the 
overall normalization) and \Ge/Gm\- The functions Hm 
and He arc replaced by the histograms obtained from 
MC simulation with the ppj selection criteria applied. 
To account for differences between the pp mass distri- 
butions of ppj events in data and MC simulation, the 
histograms Hm and He are recalculated using weighted 
events. The weights are obtained from the ratio of the pp 
mass distributions in data and simulation. In principle, 
the weights for Hm and He differ due to the different 
mass dependences of Gm and Ge- A first approxima- 
tion uses Gm = Ge- The fitted values of \Ge/Gm\ are 
then used in the next approximation to recalculate Hm 
and He- The second iteration leads to a small change 
(less than 2%) in the fitted values, and the procedure 
converges after a third iteration. 

The simulated angular distributions are corrected to 
account for the differences between data and simulation 
in particle identification, tracking, and photon efficiency. 
These corrections are discussed in detail in the next sec- 
tion. The angular dependence of the detection efficiency 
calculated with MC simulation before and after the cor- 
rections is shown in Fig. [HI The deviations from uniform 
efficiency, which do not exceed 10%, arise from the mo- 
mentum dependence of proton / antiproton particle iden- 
tification efficiency. 

The fits of the histograms to the angular distributions 



Fig.^[l + ax / {1 + bx^) , where : 



2mp GeV/c^ 



is used in the iteration procedure to calculate the weight. 
The quoted errors on \Ge/Gm\ are statistical and sys- 
tematic, respectively. The dominant contribution to the 
systematic error is due to the uncertainty in the ppn'^ 
background. 

The only previous measurement of the \Ge / Gm \ ratio 
was performed in the PS170 experiment [llj. The ratio 
was measured at five points between 1.92 GeV/c^ and 
2.04 GeV/c^ with an accuracy of 30-40% (see Fig. \TU^. 
For all points it was found to be consistent with unity. 
The average of the PS170 measurements evaluated under 
the assumption that the errors are purely statistical is 
0.90 ± 0.14. The BaBAR results are significantly larger 
for Mpp < 2.1 GeV/c^, and extend the measurements up 
to 3 GeV/c2. 



In addition, we search for an asymmetry in the pro- 
ton angular distribution. The lowest-order one-photon 
mechanism for proton-antiproton production predicts a 
symmetric angular distribution. An asymmetry arises 
from higher-order contributions, in particular, from two- 
photon exchange. Two-photon exchange is discussed 
(see, for example, Ref. [2^) as a possible source of the 
difference observed in ep scattering between the Ge/Gm 
measurements obtained with two different experimental 
techniques, namely the Rosenbluth method [24| . which 
uses the analysis of angular distributions, and the po- 
larization method |25l427l |. which is based on the mea- 
surement of the ratio of the transverse and longitudinal 
polarization of the recoil proton. 

A search for an asymmetry using previous BABAR 
e+e^ — > pp^ data [ij] is described in Ref. [2^. No asym- 
metry was observed within the statistical error of 2%. It 
should be noted that the authors of Ref. [1^ did not take 
into account the angular asymmetry of the detection cf- 
ficiency, which is seen in Fig. [9] and in a similar plot in 
Ref. [ij]. To measure the asymmetry we use the data 
with pp mass less than 3 GeV/c^. The cos Op distribution 
is fitted as described above, and the result is shown in 
Fig. [TT] Since the MC simulation uses a model with one- 
photon exchange, the asymmetry in the fitted histogram 
is due to asymmetry in the detection efficiency. To re- 
move detector effects we take the ratio of the data dis- 
tribution to the fitted simulated distribution. This ratio 
is shown in Fig. [T^ A fit of a linear function to the data 
yields a slope parameter value -0.041±0.026±0.005. The 
systematic error on the slope is estimated conservatively 
as the maximum slope given by an efficiency correction. 
The correction for the data-MC simulation difference in 
antiproton nuclear interactions (see Sec. IV1|) is found to 
yield the largest angular variation. 
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FIG. 8: The cos9p distributions for different pp mass regions: (a) 1.877-1.950 GeV/c^, (b) 1.950-2.025 GeV/c^ (c) 2.025- 
2.100 GeV/c^ (d) 2.100-2.200 GeV/c^ (e) 2.200-2.400 GeV/c^ (f) 2.400-3.000 GeV/c^ The points with error bars show the 
data distributions after background subtraction. The histograms resuh from the fits: the dashed histograms correspond to the 
magnetic form factor contributions and the dot-dashed histograms to the electric form factor contributions. 



14 



,0.22 




O 



> 



800 - 



600 - 



400 



200 



-1 



-0.5 



0.5 1 
COS 9 



FIG. 9: The angular dependence of the detection effi- 
ciency for simulated events with Mpp < 2.5 GeV/c^ before 
(open squares) and after (filled circles) correction for data- 
simulation differences in detector response. 
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FIG. 10: The measured \Ge/Gm\ mass dependence. Filled 
circles depict BABAR data. Open circles show PS170 data [ij. 
The curve is the the result of the fit described in the text. 



FIG. 11: The cos 6p distribution for the mass region from 
threshold to 3 GeV/c^. The points with error bars show 
the data distribution after background subtraction; the solid 
histogram is the fit result. The dashed and dot-dashed his- 
tograms show the contributions of the terms corresponding to 
the magnetic and electric form factors, respectively. 
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We then calculate the integral asymmetry 

a{cos9p > 0) — a{cos9p < 0) 
cr(cos 6'p > 0) -f (t(cos 6p < 0) 
= -0.025 ±0.014 ±0.003, 



(7) 



where ^(cos^p > 0) and a{cos9p < 0) are the cross sec- 
tions for e+e~ — > pp^ events vifith AIpp < 3 GeV/c^ in- 
tegrated over the angular regions vifith cos 9p > and 
cos 9p < 0, respectively. The fitted slope value and the 
integral asymmetry are consistent with zero. The value 
of the asymmetry extracted from experiment depends on 



FIG. 12: The ratio of the data distribution from Fig. [TT] to 
the fitted simulated distribution. The line shows the result of 
the fit of a linear function to the data points. 



the selection criteria used, in particular, on the effective 
energy limit for an extra photon emitted from the initial 
or final state. In our analysis, this limit is determined by 
the condition Xp < 30, and is about 100 MeV. 
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FIG. 13: The pp mass dependence of the detection efficiency 
obtained from MC simulation. 

VI. DETECTION EFFICIENCY 

The detection efficiency, which is determined using 
Monte Carlo simulation, is the ratio of true pp mass dis- 
tributions computed after and before applying the selec- 
tion criteria. Since the e^e~ — > pp^ differential cross sec- 
tion depends on two form factors, the detection efficiency 
cannot be determined in a model-independent way. We 
use a model with the \Ge/Gm\ ratio obtained from the 
fits to the experimental angular distributions (curve in 
Fig. Uni) for Mpp < 3 GeV/c2, and with \Ge/Gm\ = 1 
for higher masses. The detection efficiency obtained by 
using this model is shown in Fig. 1131 The error in the de- 
tection efficiency due to the model is determined from the 
uncertainty in the \Ge/Gm\ ratio: for Mpp < 3 GeV/c^ 
the variation of the ratio within its experimental uncer- 
tainties leads to a 1% change in the detection efficiency. 
This is taken as the model-related uncertainty. For the 
mass region above 3 GeV/c^, where the \Ge/Gm\ ra- 
tio is not measured, the model uncertainty is taken as 
the maximum difference between the detection efficien- 
cies corresponding to setting Gb = or setting Gm = 0, 
and the efficiency calculated with the model described 
above. This yields a 4% uncertainty estimate. 

The efficiency determined from MC simulation (emc) 
must be corrected to account for data-MC simulation dif- 
ferences in detector response: 

e = eMcl[{l + S^), (8) 

where the Si are efficiency corrections for each of several 
effects. These corrections are discussed in detail below 
and summarized in Table [V] 

Inaccuracies in the simulation of angular and momen- 
tum resolution and radiative corrections may account for 
data-MC differences in the fraction of events rejected 
by the requirement Xp < 30. The efficiency correction 
for this effect is estimated by comparing data and sim- 
ulated distributions for the e~^e~ — >■ fi^fi~^ process, 
which has kinematics similar to the process under study. 
An exclusive e~^e~ — > ^^/i~7 sample is selected by re- 



quiring that both charged tracks be identified as muons. 
The ratio of the number of selected muon events with 
xl > 30 and xl < 30 varies from 0.30 to 0.37 in the 
Mpp range from threshold to 4.5 GeV/c^. To charac- 
terize data-MC simulation differences in the x^ distri- 
bution, a double ratio (k) is calculated as the ratio of 
N{xl > SO)/N{xl < 30) obtained from data to the same 
quantity obtained from MC simulation. The value of the 
double ratio varies from 1.02 to 1.06 in the Mpp range 
from threshold to 4.5 GeV/c^. The efficiency correction 
Si (with i = 1) for the Xp cut is calculated as 

_ N{xl<30) + N{xl>30) 
' N{xl<SO) + KN{xl>m ' 

where N{xp < 30) and N{xp > 30) are the numbers of 
simulated ppj events with Xp < 30 and Xp > 30, re- 
spectively. The values of the efficiency correction Si for 
different pp invariant-mass values are listed in Table |Vl 

The effect of the Xk > 30 requirement is studied us- 
ing e+e^ — > J/ipl — ^ PP7 events. The J/^p yield is de- 
termined using the sideband subtraction method. The 
event losses are found to be (1.5 ± 0.4)% in data and 
(1.2 ±0.1)% in MC simulation. As the data and simu- 
lated values are in good agreement, there is no need to 
introduce any efficiency correction for the xj< > 30 re- 
quirement. The systematic uncertainty associated with 
this criterion is 0.4%. 

Another possible source of data-MC simulation differ- 
ences is is due to unreconstructed tracks. Two dominant 
effects leading to track loss in pp"f events are track over- 
lap in the DCH and nuclear interaction of protons and 
antiprotons in the material before the SVT and DCH. 

The effect of track overlap can be observed in the distri- 
bution of the parameter A(p± = ip^ — if^ , where Lp+ and 
<y9_ are the azimuthal angles at the production vertex of 
positive and negative tracks, respectively. The detection 
efficiency for simulated e^e~ — > pp^ events as a function 
of A(/3± is shown in Fig. [TH 

The z-component of the BaBAR magnetic field is in 
the direction of the positive z-axis, so that in the x — y 
plane viewed from positive z positively charged tracks 
experience clockwise bending, while negatively charged 
tracks are bent counter-clockwise. As a result, events 
with A(/J± > have a "fishtail" two-track configuration 
in which the tracks tend to overlap initially. This results 
in the dip in efficiency which is clearly seen at A(^± ^ 0.1 
rad. The ratio of the number of events with A(^± > 
to that with At/j-t < can be used to estimate the ef- 
ficiency loss due to track overlap. This efficiency loss 
reaches about 10% near pp threshold and decreases to a 
negligible level for Mpp above 2.4 GeV/c^. The effect is 
reproduced reasonably well by the MC simulation; data- 
MC simulation differences in the efficiency loss averaged 
over the mass region of maximum inefficiency, Mpp < 2.3 
GeV/c^, is about (1.2 ± 1.3)%. We introduce no correc- 
tion for this difference. For the mass region Mpp < 2.3 
GeV/c^, where the effect is large, a systematic uncer- 



16 



a 



0.3 



u 0.25 
a 
o 
■a 
o 

Q 0.2 



0.15 



0.1 



..H 



ft 



nil 



-0.5 



0.5 1 
A(|)+ (rad) 



FIG. 14: The detection efficiency for e'^e — > pp7 events as a 
function of Aip± obtained from MC simulation. 



tainty of 1.5% is assigned to the measured cross section. 

Simulation shows that nuclear interaction leads to the 
loss of approximately 6% of the e~^e~ — > ppj events. 
For data-MC simulation comparison, a specially selected 
event sample with A(A) decaying into p{p)tt is used. The 
A candidates are selected by imposing requirements on 
the pn invariant mass and A flight distance. The amount 
of material before the SVT (1.5% of a nuclear interac- 
tion length) is comparable to the amount of material be- 
tween the SVT and the DCH (1.4% of a nuclear interac- 
tion length). The probability of track losses between the 
SVT and the DCH is measured by using the A(A) sam- 
ple. The data and simulation probabilities are found to 
be in agreement for protons. A substantial difference is 
observed for antiprotons, which is consistent with a large 
(a factor of 2.6 ± 1.0) overestimation of the antiproton 
annihilation cross section in simulation. This difference 
in the antiproton annihilation cross section in data and 
simulation leads to a correction of about (1.0 ± 0.4)% to 
the detection efficiency for ppj events. 

We also incorporate a systematic uncertainty due to 
data-MC simulation differences in track reconstruction, 
which is estimated to be 0.24% per track. 

The data-MC simulation difference in particle identi- 
fication is studied using events with a. J/ip ^ pp decay. 
Due to the narrow J/tp width and hence low background, 
the number of J/ ip — >■ pp decays may be determined us- 
ing selections with cither one or two identified protons. 
The background from non- J/i/; events is subtracted using 
sidebands. The p/p identification probabilities are deter- 
mined as functions of the p/p momenta by calculating the 
ratio of the number of events with both the proton and 
the antiproton identified to the number of events with 
only one identified proton or antiproton. The ratio of 
data-MC identification probabilities is used to reweight 
selected simulated events and calculate efficiency correc- 



tions. The correction is about —(1.9 ±2.0)% and is prac- 
tically independent of pp mass. The error in the correc- 
tion is determined from the statistical uncertainty on the 
number of selected J/ip events. 

An additional correction must be applied to the pho- 
ton detection efficiency. There are two main sources of 
this correction. The first is due to data-MC simulation 
differences in the probability of photon conversion in the 
detector material before the DCH, and the second re- 
sults from the effect of dead calorimeter channels. A 
sample of e+e" — >■ events is used to determine 

the calorimeter photon inefficiency in data. Events with 
exactly two oppositely charged particle tracks identified 
as muons are selected, and a one-constraint kinematic fit 
is performed, in which the recoil mass against the muon 
pair is constrained to be zero. A tight requirement on 
the of the kinematic fit selects events with only one 
photon in the final state. The photon direction is de- 
termined from the fit, and the detection inefficiency is 
calculated as the ratio of the number of events not sat- 
isfying the E* > 3 GcV criterion, to the total number of 
selected events. The photon inefficiency obtained 

is 3.3%, to be compared to the 2% inefficiency from the 
e'^e~ — > j.i'^ii^"/ simulation. The data-MC simulation 
difference in the probability of photon conversion is also 



studied using 



H 1^ 7 events. In addition to 



two identified muons, we require that an event contain 
a converted-photon candidate, i.e., a pair of oppositely 
charged tracks with e+e~ invariant mass close to zero, 
momentum directed along the expected photon direction, 
and forming a secondary vertex well-separated from the 
interaction region. The observed data-MC difference in 
the probability of photon conversion is —(0.41 ±0.01)%. 
The data-MC diffcrcncies in the calorimeter inefficiency 
for photons and the probability of photon conversion are 
determined as functions of the photon polar angle, and 
used to reweight the simulated events and calculate effi- 
ciency corrections. The total correction due to data-MC 
simulation differences in the photon detection inefficiency 
is found to vary from of— (1.9±0.1)% near the pp thresh- 
old to — (1.7± 0.1)% at 3 GcV/c^ and higher masses. 

The quality of the simulation of trigger efficiency is 
also studied. The overlap of the samples of events satis- 
fying different trigger criteria, and the independence of 
these triggers, are used to measure trigger efficiency. The 
difference in trigger efficiency between data and MC sim- 
ulation decreases from —(0.65 ±0.20)% near pp threshold 
to -(0.13 ± 0.10)% for pp masses above 2.2 GeV/c^. An 
additional systematic uncertainty of about 0.5% is intro- 
duced to take into account the possibly imperfect sim- 
ulation of inefficiency of the offline filters that provide 
background suppression before full event reconstruction. 

All efficiency corrections are summarized in Table |Vl 
and the corrected detection efficiency values are listed in 
Table The uncertainty in detection efficiency includes 
simulation statistical error, model uncertainty, and the 
uncertainty on the efficiency correction. 
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TABLE V: The values of the difTerent efficiency corrections 
Si for pp invariant mass 1.9, 3.0, and 4.5 GeV/c^. 



effect 


5»(1.9),% 


5,(3),% 


<5,(4.5),% 


Xl < 30 


-0.5 ±0.1 


-0.9 ±0.1 


-1.5 ±0.2 


Xk >30 


0.0 ±0.4 


0.0 ±0.4 


0.0 ±0.4 


track overlap 


0.0 ±1.5 






nuclear interaction 


0.8 ±0.4 


1.1 ± 0.4 


1.0 ±0.4 


track reconstruction 


0.0 ±0.5 


0.0 ±0.5 


0.0 ±0.5 


PID 


-1.9 ±2.0 


-1.9 ± 2.0 


-1.9 ±2.0 


photon inefficiency 


-1.9±0.1 


-1.7 ±0.1 


-1.7±0.1 


trigger and filters 


-0.7 ±0.6 


-0.1 ±0.5 


-0.1 ±0.5 


total 


-4.2 ±2.6 


-3.5 ± 2.2 


-4.2 ±2.2 



THE e+e" pp CROSS SECTION AND 
THE PROTON FORM FACTOR 



VII. 



The cross section for e+e — > is calculated from the 
pp mass spectrum using the expression 



5(Mpp) 



(dN/dMpp) 

CO 

eRdL/dMp 



(10) 



where (dN/dMpp)corr is the mass spectrum corrected for 
resolution effects, dL/dMpp is the ISR differential lumi- 
nosity, s{Mpp) is the detection efficiency as a function of 
mass, and i? is a radiative correction factor accounting 
for the Born mass spectrum distortion due to emission of 
extra photons by the initial electron and positron. The 
ISR luminosity is calculated using the total integrated lu- 
minosity L and the integral over cos 9* of the probability 
density function for ISR photon emission (Eq. 



dL 

dM„i 



-((2 

TTX 



2x 



1 + B 
1 - B 



2AL 



L. 



(11) 



Q, aim determines the range of po- 
em, frame: 
20°, since 



Here B = cos6'o, and 9^ 
lar angles for the ISR photon in the e+( 
e* < e* < 180° - 61*. In our case 6* 
we determine detector efficiency using simulation with 
20° < 6*; < 160°. The values of ISR luminosity inte- 
grated over the Mpp intervals are listed in Table I VII 

The radiative correction factor R is determined from 
MC simulation at the generator level, with no detector 
simulation. The pp mass spectrum is generated using 
only the pure Born amplitude for the process e~^e~ — > 
ppj, and then using a model with higher-order radiative 
corrections included by means of the structure function 
method [l^. The radiative correction factor, evaluated 
as the ratio of the second spectrum to the first, varies 
from 1.001 at pp threshold to 1.02 at Mpp = 4.5 GeV/c^ 

The value of R depends on the requirement on the 
invariant mass of the ppj system. The value of R ob- 
tained in our case corresponds to the requirement Mpp^ > 
8 GeV/c^ imposed in the simulation. The theoretical un- 
certainty on the radiative correction calculation by the 
structure function method does not exceed 1% UM- The 
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FIG. 15: The e+e — >■ pp cross section measured in this anal- 
ysis and in other e+e~ experiments: FENI CEp , DM2[i], 
DMlH, ADONE730, BES[i, CLEO0I, NU^. The con- 
tributions of J/tp — >■ pp and il){2S) — >■ pp decays to the BABAR 
measurement have been subtracted. 
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FIG. 16: The e^e — >■ pp cross section near threshold 
measured in this analysis and in other e^e~ experiments: 
FENICE@1, DM2[|, DMl[3, ADONE730, BESQl. 



calculated radiative correction factor does not take into 
account vacuum polarization; the contribution of the lat- 
ter is included in the measured cross section. 

The resolution-corrected mass spectrum is obtained by 
unfolding the mass resolution from the measured mass 
spectrum. Using MC simulation, a migration matrix. 
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TABLE VI: The pp invariant- mass interval (Mpp), number of selected events (N) after background subtraction and mass 
migration, detection efficiency (e), ISR luminosity (L), measured cross section (cTpp), and |_Fp|, the effective form factor for 
e^e~ pp. The contributions from J/ip — >■ pp and ^(25) — >■ pp decays have been subtracted. The quoted uncertainties on A'' 
and a are statistical and systematic, respectively. For the form factor, the combined uncertainty is listed. 



Mpp (GeV/c^) 


N 




e 


L (pb-i) 


app (pb) 


\Fp\ 


1.877-1.900 


351 ± 20 ± 4 


0. 


189 ± 0.006 


2.33 


806 ± 46 ± 30 


0.424 ± 0.014 


1.900-1.925 


403 ± 22 ± 4 


0. 


178 ± 0.006 


2.52 


906 ± 50 ± 32 


0.355 ± 0.012 


1.925-1.950 


394 ± 22 ± 5 


0. 


184 ± 0.006 


2.56 


845 ± 47 ± 31 


0.309 ± 0.010 


1.950-1.975 


390 ± 22 ± 5 


0. 


186 ± 0.006 


2.60 


817 ± 46 ± 30 


0.286 ± 0.010 


1.975-2.000 


418 ± 24 ± 5 


0. 


187 ± 0.006 


2.63 


854 ± 48 ± 31 


0.281 ± 0.009 


2.000-2.025 


429 ± 24 ± 5 


0. 


192 ± 0.006 


2.67 


842 ± 48 ± 30 


0.271 ± 0.009 


2.025-2.050 


433 ± 24 ± 6 


0. 


191 ± 0.006 


2.71 


846 ± 48 ± 31 


0.266 ± 0.009 


2.050-2.075 


402 ± 24 ± 7 


0. 


197 ± 0.006 


2.75 


750 ± 45 ± 28 


0.247 ± 0.009 


2.075-2.100 


430 ± 25 ± 6 


0. 


196 ± 0.006 


2.79 


796 ± 46 ± 29 


0.252 ± 0.009 


2.100-2.125 


426 ± 25 ± 6 


0. 


195 ± 0.006 


2.83 


779 ± 45 ± 29 


0.247 ± 0.008 


2.125-2.150 


373 ± 24 ± 8 


0. 


197 ± 0.006 


2.86 


666 ± 43 ± 27 


0.227 ± 0.009 


2.150-2.175 


304 ± 22 ± 8 


0. 


192 ± 0.006 


2.90 


551 ± 41 ± 24 


0.206 ± 0.009 


2.175-2.200 


247 ± 20 ± 8 


0. 


198 ± 0.006 


2.94 


429 ± 35 ± 20 


0.182 ± 0.009 


2.200-2.225 


228 ± 20 ± 8 


0. 


198 ± 0.006 


2.98 


390 ± 33 ± 19 


0.173 ± 0.008 


2.225-2.250 


227 ± 19 ± 6 


0. 


,200 ± 0.006 


3.02 


379 ± 32 ± 16 


0.171 ± 0.008 


2.250-2.275 


139 ± 16 ± 6 


0. 


195 ± 0.006 


3.06 


234 ± 27 ± 13 


0.134 ± 0.009 


2.275-2.300 


120 ± 15 ± 6 


0. 


195 ± 0.006 


3.10 


201 ± 25 ± 12 


0.125 ± 0.009 


2.300-2.350 


173 ± 17 ± 13 


0. 


193 ± 0.005 


6.32 


143 ± 14 ± 12 


0.106 ± 0.007 


2.350-2.400 


130 ± 15 ± 13 


0. 


193 ± 0.005 


6.48 


105 ± 12 ± 11 


0.091 ± 0.007 


2.400-2.450 


143 ± 15 ± 5 


0. 


190 ± 0.005 


6.64 


115 ± 12 ± 6 


0.096 ± 0.006 


2.450-2.500 


131 ± 15 ± 5 


0. 


192 ± 0.005 


6.80 


101 ± 11 ± 5 


0.091 ± 0.006 


2.500-2.550 


111 ± 13 ± 4 


0. 


191 ± 0.005 


6.97 


84 ± 10 ± 4 


0.084 ± 0.005 


2.550-2.600 


74 ± 11 ± 4 


0. 


191 ± 0.005 


7.14 


55 ± 8 ± 3 


0.069 ± 0.006 


2.600-2.650 


55 ± 10 ± 3 


0. 


188 ± 0.005 


7.31 


40 ± 8 ± 3 


0.060 ± 0.006 


2.650-2.700 


38 ± 9 ± 3 


0. 


183 ± 0.005 


7.48 


28 ± 6 ± 3 


0.050 ± 0.006 


2.700-2.750 


50 ± 9 ± 3 


0. 


186 ± 0.005 


7.66 


36 ± 7 ± 3 


0.058 ± 0.006 


2.750-2.800 


42 ± 9 ± 3 


0. 


,184 ±0.005 


7.84 


29 ± 6 ± 3 


0.053 ± 0.006 


2.800-2.850 


25 ± 7 ± 2 


0. 


,181 ±0.005 


8.01 


18 ± 5 ± 1 


0.042 ± 0.006 


2.850-2.900 


38 ± 8 ± 2 


0. 


,174 ±0.005 


8.20 


27 ± 6 ± 2 


0.052 ± 0.006 


2.900-2.950 


28 ± 7 ± 2 


0. 


,178 ±0.005 


8.38 


19 ± 5 ± 2 


0.044 ± 0.006 


2.950-3.000 


29 ± 7 ± 2 


0. 


,170 ±0.005 


8.57 


20 ± 5 ± 2 


0.046 ± 0.006 


3.000-3.200 


25 ±12 ±9 


0. 


,168 ±0.008 


36.19 


4.2 ±2.0 ±1.6 


0.022 ± 0.007 


3.200-3.400 


36 ± 8 ± 7 


0. 


,166 ±0.008 


39.40 


5.4 ±1.2 ±1.1 


0.027 ± 0.004 


3.400-3.600 


11±4±2 


0. 


,163 ±0.008 


42.81 


1.6 ±0.6 ±0.3 


0.015 ± 0.003 


3.600-3.800 


15 ± 6 ± 2 


0. 


,167 ±0.008 


46.44 


1.9 ±0.8 ±0.3 


0.018 ± 0.004 


3.800-4.000 


1±3±2 


0. 


,168 ±0.008 


50.33 


0.2 ±0.4 ±0.2 


0.005 ± 0.005 


4.000-4.250 


4±3±2 


0. 


,164 ±0.008 


68.83 


0.3 ±0.3 ±0.2 


o.oo8t°:««t 


4.250-4.500 


3±4±2 


0. 


,160 ±0.008 


76.00 


0.3 ±0.3 ±0.2 


o.oo8«:r 



A, is obtained, which represents the probability that an 
event with true mass (Mp™*^) in mass interval j is recon- 
structed in interval i: 




The mass resolution changes from 1.5 MeV/c^ near 
threshold to 12 McV/c^ at Mpp = 3 GcV/c^ and 22 
MeV/c^ at 4.5 GeV/c^. Since the chosen mass inter- 
val width significantly exceeds the resolution for all pp 
masses, the migration matrix is nearly diagonal, with the 
values of diagonal elements ^ 0.9, and ncxt-to-diagonal 
~ 0.05. We unfold the mass spectrum by applying the in- 
verse of the migration matrix to the measured spectrum. 



The procedure changes the shape of the mass distribu- 
tion insignificantly, but increases the uncertainties (by 
w20%) and their correlations. 

After applying the migration matrix, the number of 
events in each mass interval is listed in Table IVII The 
quoted errors are statistical and systematic, respectively. 
The latter is due to the uncertainty in background sub- 
traction. The calculated cross section for e+e^ pp is 
shown in Fig.[T5]and listed in Tablc lVTl For mass intervals 
3-3.2 GeV/c^ and 3.6-3.8 GeV/c^, the nonresonant cross 
section is quoted after excluding the J/ip and tp{2S) con- 
tributions (see Sec. lVlII| . The errors quoted are statisti- 
cal and systematic. The systematic uncertainty includes 
the uncertainty on the number of signal events, detection 
efficiency, the total integrated luminosity (1%), and the 
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FIG. 17: The proton effective form factor measured in this 
analysis, in other e^e~ experiments, and in pp experiments: 
FENICEI, DMgl, DMljl, BES[|, CLEO[3, NUO, 
PS170[llfrE835[iJ E760[lJ^ (a) for the mass interval from 
pp threshold to 3.01 GeV/c^, and (b) for pp masses from 2.58 
to 4.50 GeV/c^ 



radiative corrections (1%). A comparison of this result 
with the available e'^e~ data is shown in Fig. 1151 ^-nd the 
behavior in the near-threshold region is shown in Fig. 1161 
The e~^e^ — > pp cross section is a function of two 
form factors, but due to the poor determination of the 
\Ge/Gm\ ratio, they cannot be extracted from the data 
simultaneously with reasonable accuracy. Therefore, the 
effective form factor Fp{Mpp) is introduced (Eq. (j4])), 
which is proportional to the square root of the measured 
cross section. This definition of the effective form factor 
permits comparison of our measurement with measure- 
ments from other experiments, most of which were made 
under the assumption |G£;| = \Gm\- The calculated ef- 
fective form factor is shown in Fig. [17] (linear scale) and 
Fig. [THl (logarithmic scale), while numerical values are 
listed in Table IVll These form factor values are obtained 
as averages over mass-interval width. The four mea- 
surements from PS170 [ll| with lowest mass are located 
within the first mass interval of Table I VII Consequently, 
for the mass region near threshold, where the results from 
PS 170 indicate that the form factor changes rapidly with 
mass, we calculate the cross section and effective form 
factor using a smaller mass-interval size. These results 
are listed in Table IVIIl and shown in Fig. [THl From 
Figs.[T71[IHl and[Tll it is evident that the BAR4R effective 
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FIG. 18: The proton effective form factor measured in this 
analysis, in other e^e~ experiments, and in pp experiments, 
shown on a logarithmic scale: FENlCE[il, DM2[g, DMljl, 
BESi, CLEO[3, NU[l3|, PS170[ll, E835[lll, E760[l2. 
The curve corresponds to the QCD-motivated fit described 
in the text. 



form factor results are in reasonable agreement with, and 
in general more precise than, those from previous exper- 
iments. However, in the region 1.88-2.15 GeV/c^, the 
BaBAR results are systematically above those from the 
other experiments. 

The form factor has a complex mass dependence. The 
significant increase in the form factor as the pp thresh- 
old is approached may be due to final-state interac- 
tion between the proton and antiproton p9l - l32j . The 
rapid decreases of the form factor and cross section near 
2.2 GeV/c^, 2.55 GeV/c^, and 3 GeV/c^ have not been 
discussed in the literature. The form-factor mass depen- 
dence below 3 GeV/c^ is not described satisfactorily by 
existing models (see, for example, Refs. [33l - [36j ). The 
dashed curve in Fig. [18] corresponds to a fit of the asymp- 
totic QCD dependence of the proton form factor j37l |. 



l{M^p)/M^p ^ D/{M^^\og\Myk^)l to the 
existing data with M^p > 3 GeV/c^. Here A = 0.3 GeV 
and _D is a free fit parameter. All the data above 
3 GeV/c^ except the two points from Ref. [l3| marked 
"NU" are well described by this function. Adding the 
points from Ref. [l^l changes the fit x'^/'^ from 9/16 to 
41/18, where ly is the number of degrees of freedom. The 
measurement of Ref. [l^] indicates that the form factor 
at Mpp w 4 GeV/c^ decreases more slowly than predicted 
by QCD. 
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TABLE VII: The pp invariant-mass interval (AIpp), number of selected events (A'') after background subtraction and mass 
migration, measured cross section {(Jpp), and effective form factor for e'^e^ pp {\Fp\). The quoted errors on and app are 
statistical and systematic, respectively. For the effective form factor, the combined error is listed. 



Mpp (GeV/c^) 


N 


app (pb) 


1 7"! 1 

\Fp\ 


1.8765-1.8800 


37 ± 7 ± 1 


534 ± 94 ± 39 


0.515 ±0.050 


1.8800-1.8850 


80 ± 10 ± 1 


826 ± 106 ± 42 


0.497 ± 0.034 


1.8850-1.8900 


67 ± 10 ± 1 


705 ± 105 ± 33 


0.403 ± 0.032 


1.8900-1.8950 


79 ± 11 ± 1 


886 ±121 ±41 


0.416 ± 0.030 


1.8950-1.9000 


86 ± 12 ± 1 


938 ± 128 ± 42 


0.404 ± 0.029 


1.9000-1.9050 


70 ± 11 ± 1 


785 ± 123 ± 35 


0.353 ± 0.029 


1.9050-1.9100 


80± 11 ± 1 


937 ± 135 ± 41 


0.372 ± 0.028 


1.9100-1.9150 


98 ± 13 ± 1 


1096 ± 142 ± 46 


0.390 ± 0.027 


1.9150-1.9250 


156 ± 15 ± 2 


862 ± 84 ± 32 


0.333 ±0.017 


1.9250-1.9375 


188 ± 16 ± 3 


811 ±69 ±31 


0.309 ± 0.014 


1.9375-1.9500 


208 ± 17 ± 3 


887 ± 72 ± 33 


0.311 ±0.014 


1.9500-1.9625 


181 ± 16±3 


780 ± 70 ± 30 


0.283 ±0.014 


1.9625-1.9750 


209 ± 17 ± 3 


850 ± 70 ± 32 


0.288 ± 0.013 
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FIG. 19: The proton effective form factor near pp threshold 
measured in this work and in other e'^e" and pp experiments: 
FENICE[il, DMld, PS170[ll. 



ber of J/ijj — >■ pp decays in e+e^ — > pp"f determines 
the product of the electronic width and the branching 
fraction: T{J/tp — >• e~^e~)B{J/ijj — >■ Pp)- The pp mass 
spectra for selected events in the J ftp and '0(25') mass 
regions arc shown in Fig. 1201 To determine the number 
of resonance events, each spectrum is fit with a sum of 
the probability density function (PDF) for signal plus a 
linear background term. The resonance PDF is a Breit- 
Wigner function convolved with a double- Gaussian func- 
tion describing detector resolution. The Breit-Wigner 
widths and masses for the J/i/j and ip{2S) resonances are 
fixed at their nominal values [39j . The parameters of 
the resolution function are determined from simulation. 
To account for possible differences in detector response 
between data and simulation, the simulated resolution 
function is modified by adding in quadrature an addi- 
tional (Tg to both the standard-deviation values of the 
double-Gaussian function, and introducing a shift of the 
central value of the resonance mass. The free parame- 
ters in the fit to the J/iJj mass region are the number of 
resonance events, the total number of nonresonant back- 
ground events, the slope of the background, ctg, and the 
mass shift. In the tp{2S) fit the ctg and mass shift values 
arc fixed to those obtained for the J/i/'- 



VIII. THE J/xjj AND iA(2S) DECAYS TO pp 

The differential cross section for ISR production of a 
narrow resonance (vector meson V), such as J/ip, decay- 
ing into the final state / can be calculated using [s^ 



da{s, e*) 12tt^T{V -> e+e-)B{V /) 



d cos 9X 



my s 



where my and V{y — >■ e e~) are the mass and 
tronic width of the vector meson V ^ xy = 1 — my 
and B{V — ?► /) is the branching fraction of V into the 
final state /. Therefore, the measurement of the num- 



(13) 
elec- 



The fit results are shown as the curves in Fig. [20l We 



find N 



821 ± 30 and N. 



i,{2S) 



43.5 ± 7.7: the num- 



ber of nonresonant events is 62± 11 for the 3-3.2 GeV/c^ 
mass interval, and 19 ± 6 for the 3.6-3.8 GeV/c^ inter- 
val. These values are used to extract the nonresonant 
e^e~ — > pp cross section. Since the background sub- 
traction procedure for nonresonant events (see Sec. lIVD]) 
uses events with 30 < Xp < 60, the mass spectra obtained 
with this requirement must also be fit. The numbers of 
J/ '0 and nonresonant events are found to be 40 ± 8 and 
19 ± 7. The ratio of J/ip events with 30 < Xp < 60 
to the number with Xp < 30, 0.049 ± 0.010, is in good 
agreement with the value of Ppp^ = 0.049 ± 0.001 ob- 
tained from simulation. In the '0(25') mass region, no 
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FIG. 20: The pp mass spectrum in the mass region (a) near the J/ip, and (b) near the tp{2S). The curves display the results 
of the fits described in tlie text. 



events are found with 30 < 
ing fit parameters are aa = 



Xp < 60. 
5.0 ± 1.0 



M 



The remain- 
McV/c^ and 



-(1.7 ±0.5) MeV/c2. The fitted value 



of (TG leads to an increase in simulation resolution (11 
MeV/c2) of 10%. 

The corresponding detection efficiency values are esti- 
mated from MC simulation. The event generator uses ex- 
perimental information to describe the the angular distri- 
bution of protons in J/ip and 'ip{2S) decay to pp. Specif- 
ically, each distribution is described by the dependence 



1-hacos^ Op, with a = 0.672±0.034 for J/Videcay \MbM 
and a = 0.72 ±0.13 for tp{2S) decay [4l,|4|. The model 
error in the detection efficiency due to the uncertainty of 
a is negligible. The efficiencies are found to be emc = 
0.174±0.001 for J/iP and e„c = 0.172±0.001 for i/'(2S'). 
The fractional correction for the data-MC simulation dif- 
ferences discussed in Sec. IVIl is — (3.6± 2.2)%. 

The cross section for e^e" — s- -07 — > ppj for 20° < 
61 < 160° is calculated as 



(t(20° <e*< 160° 



eRL' 



(14) 



yielding (10.4±0.4±0.3) fb and (0.55±0.10±0.02) fb for 
the J/ip and '0(25') states, respectively. The radiative- 
correction factor R = a/aBom is 1.007 ± 0.010 for the 
J/ip and 1.011 ± 0.010 for the ■0(25'), where both values 
are obtained from MC simulation at the generator level. 

The total integrated luminosity for the data sample is 
(469 ± 5) fb~^. From the measured cross sections and 
Eq. (fT3|). the following products are determined: 

r(J/Vj ^ e+e-) B{J/ij pp) = 
(11.3 ±0.4 ±0.3) eV, 
r(0(25) ^ 6+6^)6(0(25') ^ pp) = 

(0.67 ± 0.12 ± 0.02) eV. (15) 



The systematic errors include the uncertainties on the 
detection efficiencies, the integrated luminosity, and the 
radiative corrections. 

Using the nominal values for the electronic widths [s^ , 
we obtain the branching fractions 



6( J/V ^ pp) = 
6(^(25) ^ pp) 



(2.04 ±0.07 ±0.07) X 10-^ 
= (2.86 ±0.51 ±0.09) X 10""^. 



(16) 



These values are in agreement with the nominal val- 
ues dil of (2.17 ± 0.07) X 10-3 and (2.76 ± 0.12) x lO"", 
respectively, and with the recent high-precision BESIII 
result H B{J/ip pp) = (2.112±0.004±0.031) x lO^^. 



IX. SUMMARY 

The process 6"'"6- ppj has been studied for pp in- 
variant masses up to 4.5 GeV/c^. From the measured pp 
mass spectrum we extract the 6+6" — pp cross section 
and the proton effective form factor. The form factor 
has a complex mass dependence. The near-threshold en- 
hancement of the form factor observed in the PS 170 ex- 



periment jllj is confirmed in this study. There are also 
three mass regions, near 2.2 GeV/c^, 2.55 GeV/c^, and 
3 GeV/c^, that exhibit steep decreases in the form factor 
and cross section. 

By analysing the proton angular distributions we mea- 
sure the mass dependence of the ratio \Ge/Gm\ for 
Mpp from threshold to 3 GcV/c^. For masses up to 
2.1 GeV/c^, this ratio is found to be significantly greater 
than unity, in disagreement with the PS170 measure- 
ment [Tl] |. The asymmetry in the proton angular dis- 
tribution is found to be 



-0.025 ±0.014 ±0.003 
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for Mpp < 3 GeV/c2. 

From the measured event yields for e+e^ — > J/i^"/ — > 
and e~^e~ — ^ ip(2S)j — > PP7, we determine the 
branching fraction values 

B{J/ip pp) = (2.04 ± 0.07 ± 0.07) x lO'^, 
B{ip{2S) pp) = (2.86 ± 0.51 ± 0.09) x 10^1 

Our results on the cross section, form factors, and 
J/ip and ip(2S) decays agree with, and supersede, ear- 
lier BaBAR measurements [iJl . 
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